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SECTION  I 


INTRODUCTION 

This  is  the  final  report  on  Air  Force  Contract  No.  F33615- 
81-K-2003,  Project  No.  2308  for  the  experimental  and  analytical 
study  entitled,  "Laser  Velocimeter  Measurements  and  Analysis  in 
Turbulent  Flows  with  Combustion."  References  1  thru  13  are  a 
direct  result  of  the  research.  Additional  work  on  the  design, 
construction  and  preliminary  measurements  from  a  two  component 
LDV  system  are  reported  in  Section  III  of  this  report.  The 
research  .is  an  extension  of  a  previous  Air  Force  contract.  No. 
F33615-77-C-2010,  Project  No.  2308. 

The  major  results  of  this  research  are  briefly  summarized  in 
Section  II. 
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SECTION  II 


SUMMARY  OF  MAJOR  RESULTS 

1.  STUDENT  SUPPORT  AND  THESES: 

Three  MSME  theses  were  completed  with  support  from  this  con¬ 
tract.  They  are  listed  here  and  as  References  1-3. 

1.  Luchlk,  Thomas  Stephen,  "A  Laser  Velocimeter  Investigation 
of  Turbulent  Plow  in  an  Axlsymmetric  Sudden  Expansion,"  MSME 
Thesis,  Purdue  University,  May  1982. 

2.  Gould,  R.D.,  "Laser  Velocimeter  Measurements  in  an  Axisym- 
metrlc  Sudden  Expansion  With  and  Without  Combustion,”  MSME 
Thesis,  Purdue  University,  May  1983. 

3.  Durrett,  Russell  P.,  "Laser  Velocimeter  Measurements  in  an 
Axlsymmetric  Sudden  Expansion  with  a  Correction  for  Tube 
Wall  Aberrations,"  MSME  Thesis,  Purdue  University,  May  1984. 

In  addition,  the  Ph.D.  work  of  R.D.  Gould,  now  in  progress,  was 
partially  supported  from  this  contract. 

2.  REPORTS  AND  PUBLICATIONS 

In  addition  to  the  three  MSME  theses  listed  above  the  fol¬ 
lowing  ten  reports  and  publications  are  a  direct  result  of  the 
contract.  References  4,  5  and  6  are  the  interim  annual  reports 
for  the  contract  and  References  7  thru  13  are  conference  papers 
and  journal  articles  that  have  resulted  from  the  research  work. 


4.  Stevenson,  W.H.,  Thompson,  H.D.  and  Luchik,  T.S.,  "Laser 

Velocimeter  Measurements  and  Analysis  in  Turbulent  Flows 

with  Combustion,"  AFWAL-TR-82-2076 ,  Part  1,  September  1982- 

5.  Stevenson,  W.H.,  Thompson,  H.D.  and  Gould,  R.D. ,  "Laser 

Velocimeter  Measurements  amd  Analysis  in  Turbulent  Flows 

with  Combustion,"  AFWAL-TR-82-2076,  Part  II,  July  1983. 

6.  Thompson,  H.D. ,  Stevenson,  W.H.  and  Durrett,  R.P.,  "Laser 

Velocimeter  Measurements  and  Analysis  in  Turbulent  Flows 

with  Combustion,  Part  III  -  A  Correction  Lens  for  Laser 
Doppler  Measurements  in  a  Cylindrical  Tube,"  AFWAL-TR-82- 
2076,  Part  III,  July  1984. 

7.  Stevenson,  W.H.,  Thompson,  H.D.,  Gould,  R.D.  and  Craig, 
R.R.,  "Laser  Velocimeter  Measurements  in  Separated  Flow  with 
Combustion,"  Proceedings  of  the  International  Symposium  on 
Applications  of  Laser  Doppler  Anemometry  to  Fluid  Mechanics, 
Lisbon,  Portugal,  July  5-7  (1982).  Also  published  in  Laser 
Anemometry  in  Fluid  Mechanics ,  LADOAN,  Lisbon  1984  (eds.  R. 
Adrian,  D.  Durao,  F.  Durst,  H.  Mishina  and  J.  Whitelaw) . 

8.  Stevenson,  W.H.,  Thompson,  H.D.  and  Craig,  R.R.,  "Laser 
Velocimeter  Measurements  in  Highly  Turbulent  Recirculating 
F lows , "  Proceedings  of  the  Symposium  on  Engineering  Applica¬ 
tions  of  Laser  Velocimetry ,  pp.  163-170,  ASME  Winter  Annual 
Meeting,  Phoenix,  AZ,  November  1982. 

9.  Gould,  R.D. ,  Stevenson,  W.H.  and  Thompson,  ,  H.D.,  "Laser 
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Velocimeter  Measurements  in  a  Dump  Combustor,”  ASME  Paper 
83-HT-47  National  Heat  Transfer  Conference,  Seattle,  WA, 
July  1983. 

10.  Durrett,  R.P.,  Gould,  R.D. ,  Stevenson,  W.H.  and  Thompson, 
H.D.,  "A  Correction  Lens  for  Laser  Doppler  Velocimeter  Meas¬ 
urements  in  a  Cylindrical  Tube,"  AIAA  Paper  No.  84-0429, 
1984. 


11.  Stevenson,  W.H.,  Thompson,  H.D.  and  Craig,  R.R.,  "Laser 
Velocimeter  Measurements  in  Highly  Turbulent  Recirculating 
Plows,"  Journal  of  Fluids  Engineer ing  106 ,  pp-  173-180 
(1984)  . 

12.  Durrett,  R.P.,  Gould,  R.D. ,  Stevenson,  W.H.  and  Thompson, 
H.D. ,  "A  Correction  Lens  for  Laser  Doppler  Velocimeter  Meas¬ 
urements  in  a  Circular  Tube,"  AIAA  Journal  23 ,  pp.  1387-1391 
(1985) . 


13.  Durrett,  R.P.,  Stevenson,  W.H.  and  Thompson,  H.D.,  "Radial 
and  Axial  Turbulent  Flow  Measurements  with  an  LDV  in  an 
Axisymmetric  Sudden  Expansion  Air  Flow,"  to  be  published, 
ASME  Winter  Annual  Meeting,  Nov.  17-22,  1985. 

3.  MAJOR  ACCOMPLISHMENTS 

The  objective  of  this  research  has  been  to  provide  an 
increased  understanding  of  the  turbulent  flow  field  in  a  dump 
combustor.  The  progreun  has  been  primarily  experimental,  with 
some  analytical  modeling  and  comparison.  It  is  an  extension  of 


previous  work  documented  in  References  14  thru  24.  The  following 
have  been  accomplished: 

1.  A  survey  and  analysis  of  the  related  literature.  See  Refer¬ 
ences  1  thru  13. 

2.  The  identification  of  bias  errors  that  arise  in  LDV  measure¬ 
ments  and  the  development  of  measurement  techniques  to  elim¬ 
inate  those  bias  errors.  See  References  17,  20  and  23.  The 
identification  of  and  elimination  of  bias  errors  has  been  a 
matter  of  continuing  concern  throughout  both  this  and  the 
previous  research  contracts.  Major  advances  have  been  made. 
There  remain  some  experimental  situations  for  which  it  is 
neither  possible  to  apply  the  bias  elimination  procedures 
that  have  been  developed,  nor  to  accurately  compensate  for 
bias  errors  in  the  data.  Further  work  should  be  done  in 
this  area. 

3.  The  development  of  a  unique  single-component  LDV  system, 
with  variable  beam  angle  and  dual  Bragg  cell  frequency 
shifters.  The  system  could  be  rotated  around  the  beaun  axis, 
permitting  independent  velocity  component  measurements  at 
any  desired  angle  in  the  plane  perpendicular  to  the  laser 
beam  axis.  See  References  1  and  4. 

—  — — *7  — ' — '' 

4.  The  systematic  determination  of  T3,  U^,  u  »  and  u  u^^  in 

an  axisymmetric  sudden  expansion.  Measurements  were  made 
throughout  the  recirculation  region,  the  reattachment 
region,  and  in  the  recovery  region.  Logan's  method  (Ref. 
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25)  was  used  to  extract  the  desired  velocity,  turbulence  and 
shear  stress  values  from  the  individual  measurements.  Meas¬ 
ured  results  were  compared  to  computations.  See  Reference  1 
and  4. 

5.  Measurements  of  U  and  U  in  the  recirculating  and  reattach¬ 
ment  region  of  a  combusting  flow.  The  hot-flow  results  are 
compared  to  biased  and  unbiased  cold-flow  data  and  to  compu¬ 
tations.  See  References  2,  4,  7,  and  9. 

6.  The  modification  of  the  CHAMP I ON/2 /E/F IX  Computer  program  to 
handle  reacting  flow  with  the  step  geometry  used  herein. 
The  treatment  of  the  boundary  conditions  required  consider¬ 
able  modification.  An  extensive  pareunetric  study  was  con¬ 
ducted.  This  work  was  supported  in  part  by  NASA.  See 
References  5  and  26. 

7.  The  development  of  a  correction  lens  system  to  correct  for 
probe  volume  abberations  caused  by  the  circular  tube  wall. 
The  correction  lens  prevents  the  probe  volume  from  wandering 
when  the  LDV  optics  are  rotated,  and  allows  the  LDV  system 
to  be  used  to  traverse  the  flow  field  vertically.  This  in 
turn  permits  the  measurement  of  the  radial  velocity  com¬ 
ponent  and  also  allows  simultaneous  radial  and  axial  com¬ 
ponent  measurements  to  be  made  with  a  2-D  LDV  system.  See 
References  3,  6,  10,  and  12. 

8.  The  compilation  of  benchmark  cold-flow  measurements  ^ust 
downstream  of  the  sudden  expansion,  and  in  the  recirculation 
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zone.  The  correction  lens  was  used  for  these  measurements 
and  care  was  taken  to  eliminate  bias  errors.  Measurements 

were  made  to  determine  U,  u  '  and  u  u^.  See  Refer¬ 

ence  13 . 

The  development,  design,  and  construction  of  a  two-component 
LDV  system  with  improved  data  recording  and  storage  capabil¬ 
ities.  Limited  preliminary  measurements  in  a  cold-flow  rig 
have  been  made.  The  system  design  and  preliminary  data  are 
described  in  Section  III  of  this  report. 


Section  III 


THE  TWO  COMPONENT  SYSTEM 

This  section  describes  the  two-component  LDV  system  and 
presents  some  preliminary  experimental  data  obtained  in  a  cold- 
flow  sudden  expamsion  geometry. 

1.  THE  LDV  SYSTEM 

A  two-color,  two -component  LDV  system  operating  in  forward 
scatter  has  been  developed.  The  initial  application  is  simul¬ 
taneous  measurements  of  the  axial  and  radial  velocity  components 
in  an  axisymmetric  sudden  expansion  flow  with  and  without  combus¬ 
tion.  The  LDV  system  includes  Bragg  cell  modulators  in  the  four 
beam  paths  of  the  argon- Ion  laser  to  allow  a  net  frequency  shift 
of  5  MHz  in  both  the  green  and  blue  beams.  This  permits  an  unam¬ 
biguous  measurement  of  negative  velocities  and  also  eliminates 
incomplete  signal  bias.  The  green  beeun  optical  probe  volume  has 
a  waist  diameter  of  0.200  mm  and  is  approximately  2  mm  long.  The 
blue  beam  has  a  probe  volume  waist  of  0.250  mm  and  is  approxi¬ 
mately  1  mm  long.  The  scattered  light  from  the  probe  volume  is 
separated  using  a  dichroic  filter  so  that  approximately  80%  of 
each  color  passes  to  its  respective  photomultiplier  tube.  Narrow 
bandpass  filters  are  used  to  further  filter  unwanted  signals 
before  they  are  detected.  A* schematic  diagreun  of  the  LDV  system 
is  shown  in  Figure  1. 

The  flow  without  combustion  is  seeded  with  oil  ^larticles 
(DOP)  approximately  1  tim  in  diameter  generated  from  a  liquid 


atomizer  followed  by  an  evaporizat ion-condensation  unit.  The 
flow  with  combustion  will  be  seeded  with  1  ^m  diameter  aluminum 
oxide  particles  generated  from  a  TSI  model  3400  fluidized  bed 
generator  or  from  a  home-built  cyclone-  type  seeder,  depending 
upon  which  seeder  gives  the  best  results.  No  testing  has  been 
performed  using  aluminum  oxide  seed  particles  to  date.  From  some 
preliminary  testing  using  this  LDV  system  (with  seed  particles 
generated  from  the  liquid  atomizer)  maximum  data  validation  rates 
were  found  to  be  approximately  15,000  and  5,000  sounples  per 
second  for  the  green  (axial)  and  blue  (radial)  components, 
respectively.  >From  preliminary  tests,  it  appears  that  maximum 
"coincident"  data  ready  rates  range  between  1000  and  2000  per 
second.  For  this  case,  biased  velocity  data  results  as  will  be 
shown  later. 

In  an  effort  to  obtain  accurate  unbiased  velocity  data  two 
approaches  will  be  investigated.  The  first  approach  will  involve 
minor  LDV  system  modifications  with  the  goal  of  improving  the 
blue  beam  signal  quality  and  therefore  increasing  the  blue  beam 
data  validation  rate  to  >  10,000  samples  per  second.  If  this  can 
be  accomplished  the  experimental  technique  for  eliminating 
velocity  bias  by  inhibiting  the  counter  processors  for  a  fixed 
time  interval  between  samples,  approximating  equal  time  seunpling, 
will  be  used  [8,9,20,24,27,28] . 

If  the  signal  quality  of  the  blue  beam  cannot  be  improved, 
an  attempt  will  be  made  to  develop  a  velocity  bias  correction 
scheme.  The  second  approach  Involves  validating  two  proposed 


velocity  bias  correction  schemes.  The  two  correction  schemes  of 
Interest  are  the  McLaughlln-Tlederman  2-D  weighting  correction 
[29,30]  and  the  Barnett  and  Bentley  time  between  data  correction 
[31]-  The  two-component  "corrected"  data  will  then  be  compared 
to  Independent  unbiased  one-component  data  obtained  using  the 
Interval  saunpllng  technique  to  validate  and  modify  each  correc¬ 
tion  scheme. 

The  data  collection  and  processing  system  consists  of  two 
TSI  model  1990  counter-type  processors  (one  for  each  channel),  a 
TSI  model  1998  Interface  with  coincidence  timing  electronics  and 
a  PDF  11/40  mini-computer  with  DMA  capability.  A  digitized  ther¬ 
mocouple  signal  will  also  be  Interfaced  through  the  TSI  1998 
Interface  so  that  coincident  velocity-temperature  data  can  be 
obtained  (Figure  2).  With  this  system  it  Is  possible  to  acquire 
velocity  data  from  Individual  doppler  burst  signals  at  rates  up 
to  50,000  samples  per  second  (limited  by  seed  density).  The 
mini-computer  and  counter  processors  are  also  Interfaced  so  that 
sampling  can  be  controlled.  This  Is  accomplished  with  a  hardware 
clock  controlling  the  data- inhibit  hamdshaklng  signals.  The 
software  developed  for  this  experiment  Is  capable  of  burst  (meas¬ 
urement  sampled  and  stored  as  soon  as  It  is  validated)  or  Inter¬ 
val  sampling.  It  can  also  access  the  extended  memory  in  the  PDF 
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11/40  allowing  up  to  96K  words  of  data  to  be  stored  for  each  sam¬ 
ple. 


The  data  reduction  program  calculates  three  types  of  statis¬ 
tics  for  each  sample.  These  Include  the  standard  statistic,  the 
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Figure  2.  Data  Acquisition  System 


McLaughlin-Tiederni6m  2-0  weighted  statistic  and  the  time  between 
data  weighted  statistic.  Data  points  lying  outside  ±3  a  (select¬ 
able)  are  discarded  and  revised  statistics  are  calculated.  His¬ 
tograms  are  constructed  for  both  velocity  components  on  the  ter¬ 
minal  and  can  be  routed  to  the  line  printer  for  hard  copies.  The 
data  can  also  be  stored  on  the  hard  disk. 

2.  THE  TEST  RIG 


The  flow  system  utilized  in  this  experiment  is  illustrated 
in  Figure  3.  Air  is  provided  by  a  radial  fan  blower  followed  by 
a  flow  conditioning  section  consisting  of  honeycomb  flow 
straighteners .  Fuel  (gaseous  propane)  will  be  injected  in  the 
duct  immediately  following  the  blower  through  a  multi-port  mani¬ 
fold  to  give  a  homogeneous  fuel-air  mixture.  The  test  section 
consists  of  a  converging  inlet  nozzle  with  an  exit  diameter  of 
76.2  mm  followed  by  a  152.4  mm  diameter  downstream  section.  This 
inlet  was  chosen  to  give  a  uniform  inlet  velocity  profile.  The 
static  pressure  drop  across  the  nozzle  is  used  to  monitor  the 
inlet  flow  condition.  The  test  section  was  extruded  from  optical 
quality  fused  quartz  and  allows  measurements  throughout  the  flow- 
field  for  x/h  values  ranging  from  0.3  to  14.  The  test  section 
design  is  shown  in  Figure  4. 

3.  CORRECTION  LENS 

A  correction  lens  to  allow  simultaneous  axial  and  radial 
velocity  measurements  in  a  cylindrical  tube  has  been  designed  and 
fabricated  using  the  procedure  described  in  Reference  12.  The 
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plamar -concave  cylindrical  lens  that  corrects  for  the  aberration 
induced  by  the  quartz  test  section  has  a  radius  of  curvature  of 
3.3528  in,  a  thiclcness  at  its  center  of  17.93  nun  and  a  refractive 
index  of  1.52.  The  correction  lens  insures  that  the  orthogonal 
probe  volumes  (green  and  blue)  intersect  to  within  100  ^m  along 
the  length  of  the  probe  volume  and  to  within  25  fim  along  the 
diaunetex  of  the  probe  volume.  This  forces  the  scattered  light 
from  the  different  colored  probe  volumes  to  come  from  approxi¬ 
mately  the  same  point.  The  correction  lens  system  is  comprised 
of  two  lenses,  one  on  the  transmitting  side  and  one  on  the 
receiving  side  of  the  test  section  as  shown  in  Figure  5.  This  is 
due  to  the  symmetry  of  the  system.  The  lenses  have  to  be  moved 
away  from  the  test  section  as  the  measurement  point  moves  further 
off  axis.  A  ray  tracing  program  is  used  to  determine  the  place¬ 
ment  of  the  lenses  and  also  gives  the  real  probe  volume  intersec¬ 
tion  relative  to  the  beam  intersection  if  no  test  section  or  lens 
was  present.  Simultaneous  axial-radial  velocity  measurements  can 
be  made  out  to  a  non-dimensional  radius  of  approximately  85%  with 
this  system. 

4.  PRELIMINARY  EXPERIMENTAL  DATA 

A  substantial  amount  of  preliminary  two  dimensional  data 
have  been  taken  in  an  effort  to  validate  the  LDV  system,  the  flow 
characteristics  of  the  test  rig,  the  correction  lens  performance 
and  the  data  acquisition  and  reduction  software.  A  fr«e  jet  flow 
was  used  initially  and  later  replaced  by  the  more  complicated 
axlsyiranetric  sudden  expansion  flow.  As  mentioned  earlier  only 


16 


CORRECTION 

LENS 


CORRECTION 

LENS 


Figure  5.  Correction  Lens  System. 


non-reacting  flow  experiments  have  been  run  to  date.  The  follow¬ 
ing  data  set  was  acquired  using  the  2-component  LDV  and  data 
acquisition  system  described  in  the  previous  section. 

The  measurements  were  made  in  the  axisymmetric  sudden  expan¬ 
sion  flow  (Figure  4)  using  the  correction  lens  system.  The  inlet 
velocity  was  uniform  with  a  value  of  22  m/s  corresponding  to  a 

4 

Reynolds  number  of  5.5  x  10  based  on  step  height.  Seeding  was 
provided  by  the  liquid  atomizer  described  earlier.  The  burst 
sampling  technique  was  used  to  obtain  all  the  data  presented 
here.  Because  the  data  are  biased,  owing  to  the  saimpling  method 
used,  the  McLaughlin-Tiederman  2-D  correction  is  also  applied. 
This  preliminary  data  set  consists  of  measurements  across  the 
radius  of  the  test  section  at  4  non-dimensional  axial  locations, 
namely  x/H  *  0.3, 2, 4  and  6. 

Figure  6  shows  the  mean  axial  velocities  (corrected  and 
uncorrected)  at  the  four  axial  planes.  Notice  that  the  corrected 
data  has  lower  values  of  mean  velocity  in  the  shear  layer  where 
turbulence  is  high  and  has  approximately  the  same  value  as  the 
uncorrected  velocity  in  regions  of  low  turbulence.  This  is  the 
expected  result.  These  preliminary  measurements  are  presented 
only  to  demonstrate  LDV  system  performance.  They  are  known  to 
suffer  from  inaccuracy  as  a  direct  result  of  (1)  the  slowly 
oscillating  shear  layer  characteristic  of  this  flow  field  and  (2) 
the  burst  sampling  technique.  The  first  problem  occurs  because 
the  data  was  gathered  at  a  high  rate.  This  means  the  entire  sam¬ 
ple  (n  -  6000)  was  collected  in  a  very  short  time  (typically  <  2 
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Figure  7.  Normalized  Radial  Velocity 


sec)  and  the  statistical  parameters  of  the  scunpled  data  depend  on 
when  in  the  low  frequency  cycle  the  data  is  collected.  The 
observed  variation  in  the  mean  axial  velocity  due  to  this 
unsteadiness  was  up  to  10%.  This  "inaccuracy"  is  much  greater 
than  the  statistical  error  due  to  the  finite  number  of  sounples 
taken.  Typically,  in  a  well-defined  experiment,  mean  velocities 
are  repeatable  to  better  than  1%.  The  second  source  of  inaccu¬ 
racy  (the  burst  sampling  technique)  is  a  velocity  bias  problem. 
Although  the  2-D  bias  correction  was  applied,  it  was  found  to 
over  correct  the  results  in  regions  of  high  turbulence.  This  was 
confirmed  by  comparing  unbiased  one-component  (axial)  data  to  the 
2-0  corrected  data.  More  work  has  to  be  done  to  determine  the 
errors  involved  due  to  velocity  bias  and  over-correction. 

Figure  7  shows  the  mean  radial  velocity  component  at  the 
four  axial  planes.  Not  much  can  be  said  about  the  data  except 
that  the  very  low  radial  velocity  in  this  flow  field  is  respon¬ 
sible  for  more  scatter.  The  low  signal  quality  of  the  blue  beaun 
could  also  be  responsible  for  some  of  this  scatter. 

Figures  8  and  9  show  the  axial  and  radial  normalized  tur¬ 
bulence  intensities,  respectively.  The  2-D  weighted  data  show 
the  correct  trends  [9,27]  relative  to  the  uncorrected  (biased) 
data  and  reach  the  *  same  maximum  values  as  other  investigators 
have  found  [9,13].  Notice  the  shift  in  the  locations  of  peak 
turbulence  intensities  depending  on  whether  a  standard  or 
weighted  statistic  is  used.  These  figures  show  the  magnitude  of 
error  that  can  be  caused  by  biased  velocity  measurements  and 
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Figure  9 .  Normalized  Radial  Turbulence  Intens 


illustrate  that  completely  unbiased  data  are  needed  to  improve 
existing  time-averaged  turbulence  models. 

Figure  10  shows  simultaneous  Reynold's  stress  measurements 
made  with  the  coincidence  window  set  to  10  iis  (the  minimum  value 
on  the  TSI  1998  interface).  Again,  significant  differences  exist 
between  the  corrected  and  uncorrected  statistic  in  regions  of 
high  turbulence.  Also,  more  scatter  is  associated  with  this 
statistic,  probably  due  to  the  low  signal  quality  of  the  radial 
velocity  component.  However,  the  maximum  values  of  normalized 
Reynold's  stress  agree  well  with  other  experimenter's  values 
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